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Bulk and surface properties of high-quality single crystals of zirconium dodecaboride have been 
studied in the temperature range from 4.5 K up to the superconducting transition temperature 
which is found to be nearly 6.06 K. Scanning tunnelling spectroscopy data, together with dc and ac 
magnetization measurements, are consistent with the conventional s-wave pairing scenario, whereas 
they disagree in estimates of the electron-phonon coupling strength. We explain the divergence, 
supposing a great difference between the surface and bulk superconducting characteristics of the 
compound. This assertion is supported by our findings of a non-linear magnetic response to an 
amplitude-modulated alternating magnetic field, testifying to the presence of surface superconduc- 
tivity in the ZrBi2 samples at dc fields exceeding the thermodynamic critical field. 

PACS numbers: 74.70.Ad, 74.25.Op, 74.50.+r, 74.25.Ha 



Due to the unique combination of physical properties 
such as: high melting point, hardness, thermal and chem- 
ical stability metal-boron compounds have found differ- 
ent applications Q, Q • The discovery of superconduc- 
tivity in metallic MgB 2 at the unexpected T c of 39 K Q 
has caused a great interest in other transition-metal di- 
borides which show lower T c values. Related experimen- 
tal efforts have recently been extended to a wider class 
of binary boron-containing intermetallic compounds, in 
particular, to zirconium dodecaboride (HQ- 

First measurements on polycrystalline ZrBi2 (as well 
as on many other boron-rich compounds) were done by 
Matthias et al. in [|| where T c = 5.82 K was determined 
from the sharp transition in specific heat data. (Our 
estimate for T c detected by dc magnetization measure- 
ments agrees with the value of T c « 6.0 K 7]). Despite 
assertions that the ZrBi2 phase exists only at high tem- 
peratures @, large high-quality single crystals of ZrB l2 
were grown in Kiev [9j and by Leithe- Jasper et al. [lfj. 
Whereas this work reported only crystallographic 
data for the crystals, two groups 0, Q have performed 
investigations of ZrBi 2 single crystals grown in Kiev |£j 
by different physical methods. Some unexpected results 
have been derived and different conclusions relating to 
the nature of superconductivity in this compound were 
claimed. Gasparov et. al. Q observed a linear tempera- 
ture dependence of the magnetic field penetration depth 
below 3 K, which contradicts the standard BCS theory, 
and proposed a d- wave- like pairing in ZrBi2. On the 
other hand, Daghero et al. ..5J, based on resistivity vs. 
temperature data and the well-known McMillan formula 
for T c , concluded that ZrBi2 is a conventional s-wave su- 
perconductor with a zero-temperature energy-gap value 
A(0) = 0.97 meV and the ratio 2A(0)/T C = 3.64. The 
same authors Q have also measured point-contact con- 



ductance spectra that were found to be dominated by 
features typical of an s-wave superconductor but with 
the value A(T) = 0.97 meV at T nearly 4 K (the BCS 
model yields 2A(0)/T C = 4.8). The only issue that the 
two groups 0, 0| agree on, is that the single crystals 
studied are type-II superconductors with an upper crit- 
ical field H C 2(0) above 1000 Oe, as estimated in Ref. 
from resistance measurements, and in Ref. y as the field 
in which the Andreev-reflection features in conductance 
characteristics disappear. 

In this work, we have studied tunneling and mag- 
netic characteristics of ZrBi2 single crystals grown in 
Kiev. The measurements were performed for tempera- 
tures ranging from above T c down to 4.2 K, where analyt- 
ical expressions obtained within the standard Ginzburg- 
Landau approximation |ll| may be applied. In the fol- 
lowing, we address two issues: the pairing symmetry in 
ZrBi2 single crystals and whether they are really type-II 
superconductors. We have also studied surface charac- 
teristics with linear and non-linear ac magnetic response 
measurements. Chemically, dodecaborides are the ex- 
tremely stable materials in comparison with other borides 
and are characterized by a strong surface resistance to 
mechanical and chemical factors [l2j. The difference be- 
tween surface and bulk properties is discussed. 

Large rods of ZrBi 2 single crystals were grown by the 
Kiev group with typical dimensions of about 6 mm in 
diameter and up to 40 mm in length. A 10.3 x 3.2 x 1.2 
mm 3 rectangular sample was cut from the rod and was 
polished mechanically by diamond and chemically etched 
in a boiling HNO3/H2O 1 : 1 mixture for 10 min to re- 
move the Beibly layer. Room temperature X-ray diffrac- 
tion measurements performed in Kiev and in Jerusalem 
confirmed that the ZrBi2 sample is a single-phase mate- 
rial with a UB12 structure (the space group of Fm3m, 
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a = 7.407 A 0). 

The order parameter symmetry of ZrBi2 was studied 
by scanning tunneling spectroscopy measurements of cur- 
rent / - voltage V curves. In contrast to the point-contact 
technique the differential conductance dl /dV-veisus- 
V yields direct information on the local quasiparticle den- 
sity of states, and hence, on the superconductor gap sym- 
metry. The samples were carefully cleaned with ethanol 
in an ultrasonic bath just before they were mounted 
into our cryogenic home-made scanning tunneling micro- 
scope [lj]. Before inserting He exchange gas (through 
a trap) the sample space was evacuated and the device 
was dipped into a liquid-helium storage dewar. After 
sufficiently long thermalization period the sample tem- 
perature was somewhat above 4.2 K, but lower than 4.4 
K, as was indicated by a sensor placed nearby. Tunneling 
measurements were performed for junction normal-state 
resistances between 50 and 500 MOhm. I-V character- 
istics were differentiated numerically in order to obtain 
normalized spectra {dl / dV) s / {dl / dV) n , the ratios of 
differential conductances in superconducting and normal 
states. The spectra were compared with a temperature- 
smeared version of the Dynes formula that takes 
into account the effect of incoherent scattering events 
inside a superconductor, by introducing a damping pa- 
rameter r into the conventional s-wave BCS expression 
for a normalized quasiparticle density of states N(e) — 
Re[(e - iT)/y/(e-try- A(T) 2 ]. 

The well reproducible local tunneling characteristics 
(Fig. ^) did not change significantly with the tip position 
and/or the device settings. All curves demonstrated co- 
herence peaks with a pronounced minimum at V = and 
a near-gap structure symmetrical with respect to the bias 
voltage V. Gap values A(T w A.2K) were found to be 
0.97±0.01 meV (which coincides well with point-contact 
findings 5] for the same temperature) with T not exceed- 
ing 0.15 meV. A relatively small value of the damping 
factor (in particular, compared with those in Ref. |5J) is 
believed to prove the high quality of the sample surface. 
For all tip positions we have not observed any signs of a 
zero-bias peak known as a fingerprint of the d-wave pair- 
ing said, hence, may reject the assumption about an 
unusual symmetry of the order parameter in ZrBi2. 

The latter statement was also confirmed by zero-ficld- 
cooling dc magnetization M-versus-H measurements 
performed using a commercial SQUID magnetometer 
(Quantum Design MPMS) with external magnetic fields 
aligned parallel to the long axis of the sample. In the nor- 
mal state (above 6 K) the dc susceptibility \dc = M/H is 
diamagnetic as expected |7J. Fig. |2Jshows the magnetiza- 
tion curves measured at various temperatures, corrected 
by the demagnetization factor N. Usually, when the ex- 
ternal field is parallel to the long axis of a parallelepiped, 
N is small (the reduction in the H is a few percent) and 
can be neglected. However, for a sharp transition (as in 
our case) its width is of the same order as this reduc- 
tion and more detailed analysis is needed. Using the ex- 
pressions for an ellipsoid |l8| and interpolating between 
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FIG. 1: (Color online) Representative tunneling spectrum of 
ZrBi2 at T nearly 4.2 K (solid line) together with its fit to 
the temperature-smeared Dynes formula (see the text) shown 
by a dashed line, with fitting parameters A(T)=0.97 meV 
and r=0.15 meV. The inset shows the initial current-voltage 
characteristic. 
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FIG. 2: (Color online) Magnetic field dependence of the mag- 
netic moment at various temperatures: 1 - 4.5 K; 2 - 5.0 K; 
3 - 5.5 K; 4 - 5.75 K; 5 - 5.9 K. The solid line in the inset is 
an extrapolation of the thermodynamic critical field behavior 
(2) deduced from the calculated H C (T) shown by circles. 



the extreme cases that correspond to two different trans- 
verse sizes, N = 0.06 was chosen. After this correction, 
the nearly vertical magnetization curves were obtained as 
shown in Fig. [21 This vertical shape close to the critical 
field values and the tiny hysteresis loops observed (not 
shown in Fig. [2J lead to the conclusion that bulk ZrBi2 
is a type-I superconductor or, alternatively, a type-II su- 
perconductor in which the Ginzburg-Landau parameter 
k is slightly above the marginal value of k=0.71 [Til ]. 
^From the data presented in Fig. [21 we can estimate the 
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thermodynamic critical fields H Cl whose equivalent mag- 
netic energy is equal to the area under the measured M- 
versus-ff dependence J °°[-M(H)}dH = flf/(87r). Ac- 
cording to the standard BCS theory for an s-wave su- 
perconductor, the critical field curve H C (T) should satu- 
rate at low temperatures and be linear in the vicinity of 
T = T C [I3: 



flc(r) = i.735jr c (o)(i-r/r c ). 



(i) 



The latter analytical result allows us to prove the con- 
ventional symmetry of the pairing in ZrBi2, and to de- 
termine its T c value. Experimental data shown in Fig. [21 
agreed very well with the linear dependence (JQ (the 
derivative (dH c (T)/dT) T =T c was equal to -110 Oe/K) 
and an extrapolation of H C (T) to zero yielded T c = 6.06 
K. Now it is possible to use the empirical formula 



H C (T) = H C (0)[1-(T/T C 



(2) 



which interpolates the overall behavior of H C (T) between 
T = and T = T c , to estimate the zero-temperature 
thermodynamic critical field. Fig [21 (inset) yields H c (0) 
as 390 Oe. Using the density of states at the Fermi ene rgy 
calculated recently for ZrBi2 by Shein and Ivanovskii [lj| 
N(ef) = 1-687 l/(eV-cell) and the bulk energy-gap value 
of 0.97 meV from Ref. 0, one can estimate the condensa- 
tion energy within the BCS theory and to evaluate the 
thermodynamic critical field H C (Q). The value obtained 
is nearly 300 Oe which compares reasonably well with our 
estimation of H C (Q) listed above. Again, we see no rea- 
son to suppose any deviation from the conventional gap 
symmetry and the phonon origin of the Cooper pairing 
in ZrBi2. 

More problematic is the question about the electron- 
phonon coupling strength. As was pointed out by Ram- 
mer [2(j, the ratio H c (0) / (T c \dH c / dT\x=T c ) can serve as 
its indicator. In the weak-coupling limit it equals 0.58, 
as follows from Eq. (J), whereas strong-coupling effects 
lead to a significant reduction of this quantity |20|. In 
our case, this ratio equals to 0.59 which, together with 
the conclusions of the resistivity measurements 5] places 
ZrBi2 into the category of weak-coupling superconduc- 
tors. At the same time, usin g th e temperature depen- 
dence of the BCS energy gap [Tj and our value for the 
energy gap at T w 4.2 K, we obtain A(0) = (1.21 -=- 1.24) 
meV in excellent agreement with A(0) = 1.22 meV ex- 
trapolated in Ref. |5j from point-contact measurements. 
In this case the ratio 2A(0)/T C , the characteristic of the 
electron-phonon coupling strength [2l|, equals 4.75±0.10, 
indicating clearly an extremely strong electron-phonon 
interaction. According to Ref. l2ll there is a general 
trend of the 2A(0)/T C growth with increasing T c /ui\ n (ui\ n 
is a characteristic energy of lattice vibrations) and for 
2A(0)/T C =4.75, cji„ should be as great as 0.15 21]. The 
presence of such low phonon energies is doubtful, and, 
to explain this result, we assume that the order param- 
eter is higher near the sample surface. In particular, it 
should lead [2^| to a surface nucleation fields H c ^ larger 
than those expected for a uniform sample. 
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FIG. 3: (Color online) Magnetic field dependence of the dc 
susceptibility Xdc compared with the real part of the ac sus- 
ceptibility x' at the fundamental frequency Q (a) and of the 
rectified signal amplitude Act (b); T = 5.0 K. 



To study the superconducting sheath state in the field 
range H c < H < H C 3, we applied an additional small ac 
magnetic field h(t) upon the coaxial dc field H and de- 
tected the linear and nonlinear responses 23] . The mea- 
surements have been done with our original home-made 
setup 24] adapted to the MPMS SQUID magnetome- 
ter. In linear experiments the ac field h(t) = hoCOsQt 
with < ho < 0.4 Oe and H,/2tt — 1455 Hz was gener- 
ated by a copper solenoid inside the magnetometer and 
the ac susceptibility versus H was measured by the two- 
coil method. To study nonlinear characteristics of the 
sample, the ac perturbation had a form of an amplitude- 
modulated ac field h(t) = ho(l + a cos ilt) cos cot with two 
additional parameters a ~ 0.9 and w/27r = 3.2 MHz (see 
for details Ref. |24j where the same technique was applied 
to a Nb single crystal) and the amplitude of a rectified 
signal at the modulation frequency Aq was measured as 
a function of the dc field H. 

Fig. [3^, exhibits a gradual shift of the real part \' °f 
the ac susceptibility compared with \dc calculated from 
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dc ma gnet ization curves. Such behavior was explained 
in Ref. |2J as an impact of the surface superconductivity 
that appears in perfect samples in a dc field parallel to the 
sample surface and persists in a surface region adjacent to 
a vacuum interface up to a field H c ^. Nea r H C 3 is de- 
fined by a simple relation H c z — 2.38kH c The shift 
between ac and dc susceptibility curves strongly depends 
on the amplitude and frequency of the ac field h(t) and, 
because of that, any interpretation and determination of 
a critical field value H c based on ac techniques should 
be done carefully. A nonlinear response Aq was detected 
only above H c (Fig. 03) where Xdc vanishes, and hence, 
the sample bulk was in a normal state. It corresponds to 
the direct observations [^j of a nonlinear nature of the 
response wavefront when a sinusoidal ac field is applied 
to a specimen in a superconducting sheath state. Vari- 
ations of ho causes an identical behavior of the Aq in a 
Nb single crystal |24[ where the presence of surface su- 
perconductivity was proven by different authors (see, for 
example, 25] ). In our case, the ratio H c ^/H c near T c is 
about 1.8. For a vacuum interface, the Ginzburg-Landau 
parameter should exceed the marginal value of n=0.71, 
which divides type I and type II superconductors. But 
if the order parameter increases near the interface (as it 
is argued above), then the ratio of the surface nucleation 
field iJ C 3 to the thermodynamic critical field can be dra- 
matically enhanced (see Fig. 1 in Ref. |22|). Then our 
data could be interpreted with an assumption of type- 
I superconductivity. Additional experiments are needed 
to define the value of the Ginzburg-Landau parameter in 
ZrBi2 single crystals. 

In conclusion, we have proven experimentally that 



ZrBi2 single crystals are conventional s-wave supercon- 
ductors with enhanced surface characteristics. The latter 
statement can principally explain the difference in es- 
timations of the electron-phonon coupling strength and 
of critical magnetic fields between the three studies per- 
formed on the same single crystals. More investigations 
are needed to explain all findings in this non-trivial ma- 
terial, in particular, the absence of the field direction 
effect reported in Ref. H (note that surface superconduc- 
tivity should be strongly suppressed when the magnetic 
field has a component normal to the surface). We believe 
that zirconium dodecarboride is an interesting and fruit- 
ful compound for future experiments because of three 
reasons. First, due to excellent surface properties it can 
serve as a very suitable model material for studying spe- 
cific near-surface superconducting properties that are not 
yet well understood plj . Second, our experiments show 
that it is an unusual marginal superconductor near the 
border between type-I and type-II superconductors. And 
last, ZrBi2 (similar to other dodecaborides) rises above 
conventional materials due to its outstanding resistance 
to external mechanical and chemical factors. We be- 
lieve that its unique material properties and compara- 
tively simple superconducting characteristics will attract 
the attention of the applied physics community to the 
compound which can find its place among various super- 
conducting bulk applications where strong abrasion- and 
chemical resistant properties are required. 

We are thankful to R.S. Gonnelli, I.R. Shein and 
J. Kotzler for valuable discussions. The work was sup- 
ported by the INTAS program under the project No. 
2001-0617. 
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